Polarization of Prompt J/ip at the Tevatron 
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The polarization of prompt J/ip at the Fermilab Tevatron is calculated within the nonrelativistic 
QCD factorization framework. The contribution from radiative decays of P-wave charmonium states 
decreases, but does not eliminate, the transverse polarization at large transverse momentum. The 
angular distribution parameter a for leptonic decays of the J/ip is predicted to increase from near 
at pr — 5 GeV to about 0.5 at pr = 20 GeV. The prediction is consistent with measurements by 
the CDF Collaboration at intermediate values of pr, but disagrees by about 3 standard deviations 
at the largest values of px measured. 



PACS numbers: 13.85.t, 13.85.Ni, 14.40. Gx 

The production of charmonium and bottomonium 
states in high-energy collisions probes both the hard- 
scattering parton processes that create heavy quark- 
antiquark (QQ) pairs and the hadronization process that 
transforms them into color-singlet bound states. One 
particularly sensitive probe is the polarization of the 
J = 1 quarkonium states. The nonrelativistic QCD 
(NRQCD) factorization approach to inclusive quarko- 
nium production Q makes the remarkable prediction 
that in hadron collisions these states should be trans- 
versely polarized at sufficiently large transverse momen- 
tum (pt) Recent measurements at the Tevatron by 
the CDF Collaboration seem to be in dramatic contra- 
diction with this prediction || . 

As first pointed out by Cho and Wise g], the pre- 
diction of transverse polarization for 1 states at large 
Pt follows from three simple features of the dynamics of 
massless partons and heavy quarks. First, the inclusive 
production of quarkonium (or any other hadron) at suf- 
ficiently large px is dominated by fragmentation. In pp 
collisions at the Tevatron, the dominant contribution to 
the charmonium production rate at large px comes from 
gluon fragmentation ||. The gluon is almost on shell 
and thus predominantly transversely polarized. Second, 
a QQ pair with small relative momentum created by the 
virtual gluon is, at leading order in a s , in a color-octet 
3 S\ state H with the same transverse polarization as the 
gluon. Third, the spin symmetry of nonrelativistic heavy 
quarks implies the suppression of spin-flip transitions in 
the binding of the QQ pair into quarkonium. Thus, 1~~ 
states should have a large transverse polarization at suf- 
ficiently large px- A convenient measure of the polariza- 
tion is the variable a = (cr T — 2ol ) / (cry + 2ctl ) , where ax 
and <7l are the transverse and longitudinal components 
of the cross section, respectively. Beneke and Rothstein 
studied the dominant fragmentation mechanisms for o~l 
||, and concluded that, at sufficiently large px, a should 
be in the range 0.5 - 0.8. 

For charmonium production at the Tevatron, fragmen- 
tation does not yet dominate for most of the px range 
that is experimentally accessible. In order to study the 



onset of the polarization effect, it is necessary to take 
into account the fusion contributions from parton pro- 
cesses ij — > cc + k. Quantitative calculations of the po- 
larization variable a for direct ip' mesons (i.e. those that 
do not come from decays) at the Tevatron have been 
carried out by Beneke and Kramer [Q and by Leibovich 
||. They predicted that a should be small for px ^ 5 
GeV, but then should rise dramatically to 0.77 ± 0.08 
at px = 20 GeV, according to Beneke and Kramer, and 
to 0.90 ± 0.04, according to Leibovich. The CDF Col- 
laboration has measured the polarization of direct ip' || , 
but the error bars are too large to draw any definitive 
conclusions. 

The CDF Collaboration has also measured the polar- 
ization of prompt J/ip mesons |?| (i.e. those that do not 
come from the decay of B hadrons). The number of J/ip 
events is larger than for ip' by a factor of about 100, al- 
lowing a to be measured more precisely and in more px 
bins. They find that a has a positive value 0.32 ± 0.10 
in the px bin from 8 to 10 GeV. However, instead of in- 
creasing at larger px, a decreases to —0.29 ± 0.23 in the 
highest px bin from 15 to 20 GeV. Theoretical predic- 
tions of the polarization of prompt J/ip are complicated 
by the fact that the prompt signal includes J/ip mesons 
that come from decays of the higher charmonium states 
Xcij Xc2, and ip'. They account for about 15%, 15%, and 
10% of the prompt J/tp signal, respectively ||. The po- 
larization of J/ip from ip' not via Xcj is straightforward 
to calculate, since the spin is unchanged by the transi- 
tion. The polarization of J/ip from Xcj and of J/ip from 
ip' via Xcj is more complicated, because the Xcj mesons 
are produced in various spin states and they decay into 
J/ip through radiative transitions. 

In this letter, we present a quantitative analysis of the 
polarization of prompt J/ip using the NRQCD factoriza- 
tion formalism. We reanalyze the CDF data on the px 
distributions for J '/ip, Xa and "0' to determine the rele- 
vant color-octet NRQCD matrix elements (ME's). The 
cross sections for the spin states of J '/ip, Xcj, and ip' are 
calculated using these ME's and the appropriate parton 
cross sections. The cross sections for the spin states of 
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the Xcj required the calculation of new parton cross sec- 
tions, which will be published elsewhere pC] . The vari- 
able a is then obtained by combining these cross sections 
with the appropriate branching ratios into longitudinally 
polarized J/ip {ipL,)- 

The NRQCD factorization formula for the differential 
cross section for the inclusive production of a charmo- 
nium state H of momentum P and spin quantum number 
A has the schematic form 

da H ^ = da c5 ^ (0^ (P) >, (1) 

where the summation index n runs over all the color 
and angular momentum states of the cc pair. The cc 
cross sections da CCn can be calculated using perturba- 
tive QCD. All dependence on the state H is contained 
within the nonperturbative ME's (0^ A ' P '). In general, 
they are Lorentz tensors that depend on the momentum 
P and the polarization tensor of H\. The Lorentz in- 
dices are contracted with those of da CCn to give a scalar 
cross section. The symmetries of NRQCD can be used 
to reduce the tensor ME's (O r f a(P) ) to scalar ME's (0%) 
that are independent of P and A. Thus one may cal- 
culate the cross section for polarized quarkonium once 
the relevant scalar ME's are known. A nonperturbative 
analysis of NRQCD reveals how the various ME's scale 
with the typical relative velocity of the heavy quarks. 
It also gives exact and approximate symmetry relations 
that can be used to simplify the ME's. The most im- 
portant ME's for the production of J/ip — ipi)-S) or 
ip' = ip(2S) can be reduced to one color-singlet pa- 
rameter (O^ ( 3 /Si)) and three color-octet parameters 

(<# (bS) ( s £i)>, (0f" S) (%)), and (of nS \ 3 P )). The 
most important ME's for Xcj production can be reduced 
to a color-singlet parameter (O* c0 ( 3 -Po)) and a single 
color-octet parameter (Og c0 ( 3 Si)). The ME's enumer- 
ated above should be sufficient for a calculation of the 
polarization of prompt J/ip. 

In pp collisions, the parton processes that dominate 
the cc cross section depend on px- If pr is of order m c , 
those which dominate are fusion processes, whose contri- 
butions can be expressed as 

<: {p) - h/v ® fj/P ® d °7 (p) (°n a(p) >> (2) 

where fi/ p {x,p) and fj/p{x,n) are parton distribution 
functions (PDF's) and a sum over the partons i,j is im- 
plied. The leading-order parton cross sections da are 
proportional to a 3 (p). These cross sections are given in 
Refs. Q and JllJ for all the relevant cc spin states with 
the exception of color-singlet 3 Pj states, which required 
a new calculation. For pt 3> m c the parton cross sec- 
tions are dominated by fragmentation processes with the 
scaling behavior da/dp T ~ l/py- These contributions 
can be expressed as 

< A(P) = f i/p ® Sm ® **T /Z) ® D t n <°« A(P) )' ( 3 ) 



where D^ Cn (z , pf Y ) is a fragmentation function (FF). We 
use a common renormalization and factorization scale 
|U for fi/p, fj/p, and da, but we allow /if r to be differ- 
ent. The momentum k of the fragmenting parton is de- 
noted by P/z in (||). However, it is inconsistent to set 
fan — p^/ z> since the parton is massless while P 2 = 4m 2 . 
We choose fc M so that z is the fraction of the light-cone 
momentum of the parton k that is carried by the cc pair 
in the parton CM frame. The covariant expression is 
k^ = [(A + K ■ P)P^ - P 2 K^\ /(2zA), where is the 
total momentum of the colliding partons i and j, and 
A=[(K-P) 2 -K 2 P a ]V 2 . 

In order to predict the polarization of prompt J / ip at 
the Tevatron, we need values for the scalar ME's. The 
color-singlet ME's (Of {nS) ( 3 S*i)) and (0^"( 3 P a )) can be 
determined phenomenologically from the decay rates for 
tp(nS) — ► £ + £~ and \c2 — > 77 |0- Using the vac- 
uum saturation approximation and spin symmetry in the 
NRQCD factorization formulae and including NLO QCD 
radiative corrections [ fl3|| , we obtain the values in Table 
[|. The errors come from the experimental errors in the 
decay rates only. 

The color-octet ME's are phenomenological parame- 
ters that must be determined from production data. To 
predict the polarization at the Tevatron, it is preferable 
to use ME's extracted directly from Tevatron data in or- 
der to cancel theoretical errors associated with soft gluon 
radiation. There have been several previous extractions 

of the color-octet ME's f||l|-|J] from the CDF data on 
the pt distributions of J/ip, Xo and ip' ||. We carry out 
an updated analysis largely following the strategy used 
in Ref. p6[ . In the fusion cross section (g), we include 
the parton processes ij — > cc+k, with i,j = g,q,q and 
q = u,d,s. In the fragmentation cross section (g), we 
include only the g — ► ccs( 3 S'i) term, since this is the only 
fragmentation process for which D^ n is of order a s . The 

FF Df 8( Sl) is evolved in /Zf r using the standard homo- 
geneous timelike evolution equation. The effects of the 
violation of the phase-space constraint /if r > Am 2 / z are 
negligible at the Tevatron due to the rapid fall-off of the 
Pt distribution [ fL7| . 

We consider two choices for the PDF's: MRST98LO 
as our default and CTEQ5L for comparison We 
evaluate a s from the one-loop formula using the value 
of Aqcd appropriate for the PDF set [Tq| . We set 
= (4m 2 + Pt) 1 ^ 2 an d m c = 1.5 GeV. The cross 
section for tp(nS) depends on the linear combination 
M r = (O 8 ( 1 5 )) + r(O$( 3 P )}/m 2 , where r varies from 
about 3.6 at pt = 5.5 GeV to about 3.0 at pt = 18 GeV, 
so we can only determine M r at some optimal value of r. 
We determined (0 8 ( 3 Si)) and M r for ip(nS) by fitting 
the pt distributions from CDF following the strategy in 
@. We determined (0% c0 ( 3 Si)) by fitting the p T distri- 
bution for Xc together with the constraint from the pre- 
liminary CDF measurement of a Xrl /a Xc2 M . Our values 



2 



for the color-octet ME's are summarized in Table ffl. The 
error bars take into account the statistical errors only. 
Our default ip' color-octet ME's agree within errors with 
those of Ref. used by Leibovich Q and with those for 
2 of the 3 PDF sets used by Bencke and Kramer . Our 
default J/ip color-octet ME's agree within errors with 
those of Ref. @. Our value for (Og /V '( 3 S'i)) is about a 
factor of 3 smaller than in Ref. JiJ , while Mr is about 
a factor of 2 larger. 

We can calculate the cross sections for the polarized 
states H\ using the scalar ME's in Table |[ The cross 
section (Q) can be reduced to an expression linear in the 
scalar ME's, with coefficients that involve the polariza- 
tion tensor of H\. In the channel ccs( 3 Si) — » ip\(nS), 
we interpolate between the fusion cross section at low pt 
and the fragmentation cross section at high pt using the 
prescription 

da H > = da? u * x (daf/ [/x fr = fi]/da^ [ Mfr = 2m c ]) . (4) 

We proceed to summarize our calculation of the errors 
in ctl and <tt- The errors in the ME's in Table Q are 
taken into account. We take the central values of \i and 
m c to be /it = (4m;? + Pt) 1 / 2 and 1.5 GeV and allow 
them to vary within the ranges \[It ~ 2/ir and 1.45 - 
1.55 GeV, respectively. We take MRST98LO as our de- 
fault PDF, and we treat the difference between it and 
CTEQ5L as an error. The cross section gl for ip(nS) is 
sensitive to a different linear combination of (Os^Sq)} 
and (Os( 3 Po)) than appears in M r . We take this into 
account by expressing the cross section as a function of 
M r and x — (OsC Sq)) /M r , taking the central value of x 
to be \i and allowing x to vary between and 1. 

We first consider the polarization of direct ip' , since 
it is not complicated by feeddown from higher charmo- 
nium states. The polarization variable a measured by 
the CDF Collaboration j|] describes the angular distri- 
bution of leptons from the decay of the ip' with respect 
to the ip' momentum in the hadron CM frame. The co- 
variant expression for the polarization vector of ip' L is 
(P 2 Q^ -P-QP> 1 )/{\/P 2 A), where Q = p + p is the total 
hadron momentum and A = [(P ■ Q) 2 — P 2 Q 2 ] 1 / 2 . In 
Fig. 1(a), we compare our result for a as a function of 
Pt with the CDF data || and with previous predictions 
from Refs. @ and || . We present our result in the form 
of an error band obtained by combining in quadrature 
all the errors described above. The most important er- 
rors are those from the ip' ME's, the PDF's, and x. The 
error bars in the CDF data are too large to draw any 
definitive conclusions. Our result for a is close to the 
prediction of Leibovich ]8| , and significantly larger than 
that of Beneke and Kramer JtJ . Their calculations differ 
in the treatment of terms of order a 2 in the gluon FF 
. Beneke and Kramer included these terms in <tl but 
neglected them in <jt, while Leibovich neglected them 
in both ot and ot- We have adopted the strategy of 



Beneke and Kramer, since these terms give a significant 
increase in ol at large pt but have only a small effect on 
ot- Although this tends to decrease a, our smaller value 
of M r tends to increase a, and the net result is close to 
the prediction of Leibovich. 

We next consider the polarization variable a for 
prompt J /ip. The prompt cross section <jt + o\l is the 
sum of the direct cross section for J/ip and the cross sec- 
tions for Xcj and ip' weighted by the branching fractions 
Bh^.j/iP- The prompt longitudinal cross section at, is 
the sum of the direct cross section for ip^ and the cross 
sections for each of the spin states Xc.j(X) and ip' x weighted 
by B H _,j/^ and by the probability Ph x ^ l for the po- 
larized state to decay into ipL- The observed transitions 
of ip' to J/ip involve no spin flips, so that P^i ^,^ L is 1 
for ip' and for ip' ±1 . For the radiative decay of Xcj(A) 
into J/ip, the probability P XcJ(x) ^ L is ^ for Xcth \ f° r 

Xci(±i), | for Xc2(o), \ for Xc2(±i), and for the other 
spin states [Q. In Fig. 1(b), we compare our result for 
a as a function of pt with the CDF data || . The shaded 
area indicates the error band obtained by adding the er- 
rors in quadrature. The most important errors are those 
from the PDF's, the J/ip ME's, and x. Our result for 
a is small around pt = 5 GeV, but it increases with pt 
to a value around 0.5 at pt — 20 GeV. Our result is in 
good agreement with the CDF measurement at interme- 
diate values of pt, but it disagrees by about 3 standard 
deviations in the highest pt bin. The three solid lines 
in Fig. 1(b) are the central curves of a for direct J/ip, 
J/ip from Xcj, and J/ip from ip'. The a for direct J/ip is 
smaller than that for direct ip' , because (Os( 3 Si)) is com- 
parable for J/ip and ip' , while M r is significantly larger for 
J /ip. In the moderate-pT region, the contributions from 
ip' and from Xc add to give an increase in the transverse 
polarization of prompt J/ip compared to direct J ftp. In 
the high-pT region, the contributions from ip' and \c tend 
to cancel. The prediction of Beneke and Kramer for a 
for direct J/ip is identical to their prediction for direct 
ip' in Fig. 1(a). At pt — 20 GeV, it is significantly larger 
than our prediction for direct J/ip. The difference comes 
from our smaller value for (Og^( 3 S'i)) and our larger 

value for Af/^. Beneke and Kramer's prediction for a 
for J/ip from ip' would be significantly lower than our re- 
sult in Fig. 1(b), but it would have a small effect on the 
value of a for prompt J/ip. The discrepancies between 
their predictions and ours could be eliminated by more 
accurate data on the J/ip and ip' cross sections, which 
would decrease some of the ambiguities in the analysis. 

The CDF measurement of the polarization of prompt 
J/ip presents a serious challenge to the NRQCD factor- 
ization formalism for inclusive quarkonium production. 
There are many effects that could change our quantita- 
tive prediction for a, such as next-to- leading order ra- 
diative corrections, but the qualitative prediction that a 
should increase at large pt seems inescapable. In Run II 
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PDF {Of+fSi)) (0 8 J/ *( 3 Si)> M 3 J/ / (0f( 3 Si)> (Ot'( 3 S 1 )) Mf 5 (O*°°( 3 P )) (O^Sr)) 

MRST98LO 1.3 ±0.1 4.4 ±0.7 8.7 ±0.9 6.5 ± 0.6 4.2 ±1.0 1.3 ± 0.5 8.9 ±1.3 2.3 ± 0.3 

CTEQ5L 1.4 ±0.1 3.9 ± 0.7 6.6 ± 0.7 6.7 ±0.7 3.7 ±0.9 0.78 ± 0.36 9.1 ± 1.3 1.9 ±0.2 

unit GeV 3 10 _3 GeV 3 10 _2 GeV 3 10 _1 GeV 3 10 _3 GeV 3 10 _2 GeV 3 10 _2 GeV 5 10 _3 GeV 3 



TABLE I. NRQCD matrix elements. The error bars take into account the statistical errors only. 




of the Tevatron, the data sample for J/tp should be more 
than one order of magnitude larger than in Run I, allow- 
ing the polarization to be measured with higher precision 
and out to larger values of px- If the result continues to 
disagree with the predictions of the NRQCD factoriza- 
tion approach, it would indicate a serious flaw in our 
understanding of inclusive charmonium production. The 
predictions of low-order perturbative QCD for the spin- 
dependence of cc cross sections could be wrong, or the use 
of NRQCD to understand the systematics of the forma- 
tion of charmonium from the cc pair could be flawed, or 
m c could simply be too small to apply the factorization 
approach to the charmonium system. 
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